Controlling a diatomic shape resonance with non-resonant light 
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| | A (diatomic) shape resonance is a metastable state of a pair of colliding atoms quasi-bound by the centrifugal barrier imposed 
.by the angular momentum involved in the collision. The temporary trapping of the atoms' scattering wavefunction corresponds 
^| to an enhanced atom pair density at low interatomic separations. This leads to larger overlap of the wavefunctions involved in a 
s^.molecule formation process such as photoassociation, rendering the process more efficient. However, for an ensemble of atoms, 
jthe atom pair density will only be enhanced if the energy of the resonance comes close to the temperature of the atomic ensemble. 
^Tj .Herein we explore the possibility of controlling the energy of a shape resonance by shifting it toward the temperature of atoms 
'confined in a trap. The shifts are imparted by the interaction of non-resonant light with the anisotropic polarizability of the atom 
pair, which affects both the centrifugal barrier and the pair's rotational and vibrational levels. We find that at laser intensities of 

up to 5 x 10 9 W/cm 2 the pair density is increased by one order of magnitude for 87 Rb atoms at 100 fiK and by two orders of 

C| magnitude for 88 Sr atoms at 20 jiK. 
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1 Introduction 
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The quest for translationally ultracold (T < 100 ;llK) 
molecules in their absolute internal ground state has been re- 
warded recently when ultracold atoms were associated via 
a magnetic-field controlled Feshbach resonance followed by 
Stimulated Raman Adiabatic Passage.^ 2 - Ultracold molecules 
(— | were also produced by associating atoms using laser light— and 
brought into their ro vibrational ground state by spontaneous 

1 1 or stimulated emission.ir£ Feshbach or magneto-association 

I requires non-zero nuclear spin - a hyperfine manifold - of 
^> the colliding atoms.— It is efficient when the collision energy 
i— I of the atoms is of the order of the hyperfine splitting. This 
corresponds to nano-Kelvin temperatures which are reached 
via evaporative or sympathetic cooling, leaving typically only 
about 10 4 — 10 s atoms that can be associated. Photoassocia- 



o 



in 
o 



- 1—1 

x 



tion^a relies only on the presence of optical transitions which 
are usually abundant, and is not tied to any particular temper- 
ature regime (although the specific photoassociation mecha- 
nisms may differ at high~ and very low— temperatures). In 
the ultracold domain, photoassociation has been often imple- 
mented in combination with magneto-optical traps (MOTs), 8 
which hold up to 10 10 atoms at temperatures ranging between 
.'10— 100 jxK. However, only a small fraction of the MOT 
atoms are at internuclear separations, R, amenable to photoas- 
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sociation, i.e., the density of atom pairs within the required 
range of R is tiny. This reflects the fact that photoassocia- 
tion cannot achieve the phase space compression required to 
transform far-apart atoms into tightly or even weakly bound 
molecules. As a result, the number of ground state molecules 
produced via photoassociation has been limited to only about 

1Q 4 5,6 

A remedy for the limited efficiency of photoassociation 
could be found within the field of coherent control, which 
employs quantum interference to constructively enhance a de- 
sired outcome of a process while destructively suppressing 
all its undesirable alternatives. 1Q i n Although coherent control 
has been successful for unimolecular processes such as pho- 
toionization or photodissociation, 12 ' 13 controlling a binary re- 
action has remained an open challenge. This has been mainly 
due to the fact that the initial state of the reaction consists of 
an incoherent mixture of scattering states with random rela- 
tive phases. However, resonances can, in principle, endow the 
mixture of the initial scattering states with a single phase.— 

The way toward solving the problem of limited photoasso- 
ciation efficiency may be paved by considering photoassoci- 
ation as a coherent control problem. Coherent control em- 
ploys quantum interference to constructively enhance a de- 
sired outcome of a process while destructively suppressing 
all its undesirable alternatives. 1Q i u Although coherent control 
has been highly successful for unimolecular processes such as 
photoionization or photodissociation, 12,13 controlling a binary 
reaction such as photoassociation remains an open challenge. 
This is mainly due to the fact that the initial state of the reac- 
tion consists of an incoherent mixture of scattering states with 
random relative phases. However, resonances can, in princi- 
ple, endow the mixture of the initial scattering states with a 
single phase.— 
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The importance of resonances increases in the ultracold 
regime where, along with tunneling, they dominate the quan- 
tum dynamics. The presence of a Feshbach resonance has 
been predicted to significantly enhance the photoassociation 
yield in Feshbach optimized photoassociation (FOPA).— Sim- 
ilarly, an enhanced microwave absorption is expected near a 
Feshbach resonance.— However, Feshbach enhancement of 
photoassociation is restricted to atoms with non-zero nuclear 
spin and is the more efficient the lower the temperature. Al- 
though electric field-induced resonances^"— should also ic- 
nrease the photoassociation efficiency, the electric fields re- 
quired to achieve a significant enhancement are currently ex- 
perimentally unfeasible. Shape resonances that occur when a 
scattering state becomes trapped behind the centrifugal bar- 
rier for partial waves with / > were found to yield enhanced 
photoassociation rates.^2r— However, due to the rotational ex- 
citation involved in generating the barrier, the lowest energies 
at which shape resonances occur correspond typically to tem- 
peratures of a few milli-Kelvin. Therefore, the thermal weight 
of a shape resonance in a much colder MOT or optical trap is 
quite small. 

Herein, we study the possibility of enhancing the thermal 
weight of a scattering resonance by shifting its position (with 
respect to the trap energy) by applying a non-resonant radia- 
tive field. The interaction of such a field with the anisotropic 
polarizability of the atom pair is of a universal character, in- 
dependent of any particular energy level structure, frequency 
of the light (as long as it is non-resonant), or the presence of a 
permanent dipole moment. Non-resonant or far-off-resonant 
light was used to manipulate molecular alignment—"— and 
has been predicted to shift rotational and vibrational levels — 
or to modify intermolecular interactions 2 - 7 , It has also been 
predicted to cause peculiar effects in atomic Bose-Einstein 
condensates such as gravitational self -binding and supersolid- 
like structures.—"— The nonresonant polarizability interac- 
tion creates an effective rotational barrier for weakly bound 
molecules; and the subsequent 'shaking' of the molecules im- 
parted by the nonresonant light could be employed to recover 
the vibrational probability density distribution of the corre- 
sponding bound atom pair. 31 > 32 The intuitive picture of an ef- 
fective centrifugal barrier provided an impetus for the present 
study: if the non-resonant light creates (for J = 0) or modi- 
fies (for J > 0) the centrifugal barrier, it will also affect the 
position of a shape resonance. Application of a non-resonsant 
field should therefore allow to control a shape resonance. The 
core of our present study extends the previous work on shift- 
ing weakly bound diatomic levels to shifting of quasi-bound 
scattering states. 

The ubiquitous diatomic low-/ shape resonances are 
straightforward to predict from the basic scattering properties 
of ultracold atoms. 22 i 33 i 34 Specifically, we consider a of-wave 
shape resonance in 87-rubidium 2 • and a g-wave shape reso- 



nance in 88-strontium. Rubidium is the drosophila of ultra- 
cold physics, whose potential energy curves and polarizability 
anisotropy are accurately known. 88 Si"2 features in proposals 
to test the time dependence of the electron-to-proton mass ra- 
tio.— With no Feshbach resonances present, 88 Sr2 has been 
produced by photoassociation. 36 > 37 We compare the prospects 
for controlling the shape resonances of rubidium and stron- 
tium with non-resonant light and draw general conclusions 
about the applicability and efficiency of this control scheme 
for producing diatomics via photoassociation. 

2 Theoretical framework 

2.1 Description of a thermal cloud of atoms 

The initial state of an ensemble of atoms held at thermal equi- 
librium in a trap of temperature T is described by the canonical 
density operator, 
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Equation (Q]i disregards the effects of quantum statistics, 
which are negligible at typical MOT temperatures. In a dilute 
gas where three-body and higher order interactions can be ne- 
glected, the Hamiltonian, Hjy, consists only of single-particle 
(T, Vtrap) an d two-body Vy(Ry) operators. The /V-particle 
density operator, Pt.n, is then given in terms of 7V 2 -times the 
pair density operator, pj- 2 .— Assuming that the center of mass 
and internuclear degrees of freedom can be separated - as is 
the case, e.g., in a harmonic trap - the center of mass motion 
can be integrated analytically. By representing the Hamilto- 
nian for the internuclear motion on a finite-size coordinate- 
space grid with variables (R, 9,(p), and making use of the az- 
imuthal symmetry, the density matrix of the initial state can 
be constructed in terms of the eigenfunctions, (p n j(R,0), of 
the pair Hamiltonian, H2, with eigenvalues E„j,— 
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The grid needs to be sufficiently large in R to approximate 
the scattering continuum well.— A small number of partial 
waves is sufficient to ensure convergence of the short-range 
part of pj 2 (R, 0), since only a few partial waves are thermally 
populated at very low temperatures. 38 

Formally, the time evolution of the cloud of atoms is deter- 
mined by the dynamics of the pair density operator, 



p T2 (t)= e T> Hl ' p T ^(t = 0)e* 



(3) 



under the assumption that on the timescale of t, the time evo- 
lution is unitary, i.e., no dissipative mechanisms are present. 
Any thermal expectation value is obtained as (A(f)) = 
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Tr[A/$r,2(f)]- F° r exam pl e i the photoassociation probability 
is given by (P e (tf)), where P e denotes the projector onto the 
electronic state e which is populated by photoassociation and 
tf is some final time. The unitary time evolution in Eq. (O 
implies that it is not necessary to solve the Liouville-von Neu- 
mann equation for the density operator explicitly. Rather, a 
separate propagation for each eigenstate in Eq. (ff) is sufficient 
to calculate thermal expectation values.— In particular, if only 
the short range part of the initial thermal density is probed, 
only few partial waves suffice for convergence. This is nu- 
merically much more efficient than solving the Liouville-von 
Neumann equation since the number of terms in the sum of 
Eq. (|2} is relatively small due to the narrow thermal width of 
an initial state at low temperature. 

2.2 Interaction of a diatom with non-resonant laser light 

The Hamiltonian of an atom pair in its electronic ground state 
in the presence of a non-resonant laser field can be written as 

Hi = T R + ^ + V g (R) - — (Aa(R)cos 2 + a ± (R)J , 

(4) 

where the first and second terms denote the vibrational and 
rotational kinetic energies, respectively, and V^(R) the field- 
free ground electronic potential energy curve. The last term 
of Eq. represents the interaction with a non-resonant laser 
field, with / the laser intensity and Aa(R) = a,\\ (R) - a ± (R) 
the polarizability anisotropy, given in terms of the perpendic- 
ular and parallel molecular polarizability components, a±(R) 
and OCjj(R). Forhomonucleardimers, the behaviour at large in- 
ternuclear distances is given by Silberstein's expansion,^^ 3 - 

<Xi(R) « 2ao-2c^/R 3 +2c^/R 6 , 

«H(R) » 2o + 4a 2 /R 3 + 8a 3 /R 6 , (5) 

where ao is the atomic polarizability. Hamiltonian © is de- 
rived by assuming the frequency of the laser to be far from 
any resonance and larger than the inverse of both the pulse 
duration and the rotational period. In this case, a two-photon 
rotating-wave approximation averaging over the rapid oscilla- 
tions of the non-resonant field can be applied. — Furthermore, 
far from resonances, the frequency-dependent molecular po- 
larizability approaches its static value, which allows to cast 
Eq. © in the static polarizability limit. The second-order na- 
ture of the light-matter interaction is reflected by the intensity 
(not the field amplitude) and cos 2 operator occurring in the 
last term of Eq. ©. Thus the energy of the field-dressed eigen- 
states always decreases with increasing field intensity, i.e., the 
states are high-field seeking. 

Since an external field defines a preferred direction in space, 
the symmetry of the corresponding Hamiltonian is reduced 
compared to the field-free case. In the absence of a field, 



a bound molecular state is characterized by its vibrational, 
rotational, and magnetic quantum numbers (v,J,M). This 
carries over to unbound box-discretized states approximat- 
ing the scattering continuum above the dissociation limit, i.e., 
(n,J,M) where each translational "quantum number" repre- 
sents a range of scattering states with collisional energies close 
to E n .— Due to the azimuthal symmetry about the laser po- 
larization axis, the light-matter interaction depends only on 
the polar angle 6. As a consequence, by turning on the non- 
resonant laser field, hybridization of the rotational motion 
takes place and only the magnetic quantum number M remains 
conserved. The field-free degeneracy of states with the same 
rotational quantum number J but different M is lifted due to 
the interaction with the laser field. Herein, we focus on states 
with M = 0, since for a given J the effect of the non-resonant 
field increases with decreasing M and is largest for M = 0. 
Thus, the sum in the density matrix (O describing the initial 
state is restricted to states with M = 0, setting the degeneracy 
factor (27 + 1) to 1. 

Discrete basis set methods are employed to represent the 
vibrational and rotational degrees of freedom of Hamiltonian 
©. For the radial part, a mapped Fourier grid is used with the 
grid step set to be proportional to the local de Broglie wave- 
length. 45-47 This allows us to use large enough discretization 
boxes to properly describe the part of the scattering continuum 
that is relevant at ultracold temperatures. The rotational de- 
gree of freedom is treated by a discrete variable representation 
in terms of Legendre polynomials, taking into account that M 
is conserved. 48 ' 49 At large intensities, more partial waves than 
those that are thermally populated initially will come into play, 
cf. Eq. (|2]). Convergence of our calculations with respect to 
the number of Legendre polynomials is ensured for the largest 
intensity employed (Jj najc « 14 compared with jf nax = 3 or 4 in 
the field-free case). 

As evident from Eq. ©, the interaction with the non- 
resonant laser light couples R and 9, i.e., it affects both vi- 
brational and rotational motion of the atom pair. In order 
to ease interpretation, it is expedient to disentangle the two 
effects by making use of effective one-dimensional models. 
The effect on the rotational motion is captured by an effective 
quantum number, J*, determining the rotational barrier in an 
effective one-dimensional (vibrational) model.~ If the non- 
resonant laser field is kept on, the effect of the field on the 
angular motion at a given internuclear separation needs to be 
accounted for by including the admixture of different partial 
waves and of the alignment. We therefore generalize previous 
treatments 23 i 24 i 26 i 32 i 50 i 51 to the latter case. In order to obtain 
an effective one-dimensional model in the presence of the non- 
resonant light, we first solve the angular part of the eigenvalue 
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problem for a fixed R, 



|2 27tl / 

[ Aa(R) cos 2 9 + a ± (R) 
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3>j(9;R,I) 
= Ej(R,I)<t>j(9;R,I). (6) 



The eigenfunctions, <Pj(9;R,I), and eigenvalues, Ej(R,I), of 
this equation depend parametrically on the radial coordinate 
R and on the laser intensity /. The index j is a label which 
is related to the field-free rotational quantum number. The 
full wavefunction can be expanded in terms of these angular 
wavefunctions, 



(7) 



Inserting this expansion into the time-independent 
Schrodinger equation for Hamiltonian ©, multiplying 
by <&*j(0;R), integrating over the angular coordinate 6, 
and neglecting all non-adiabatic coupling matrix elements, 
Jo<Z>*(9;R,I)t R <i> j <(6;R,I)smGde, the following effective 
one-dimensional Schrodinger equation for the vibrational part 
is obtained, 

[t R + V g (tL)+Ej(tL;I)] y nj {R;I)=E n y nj {R;I). (8) 

Here, V g (R) +Ej(R;I) is the effective potential including the 
effect of non-resonant light. Note that although the index « 
is used as a label indicating a scattering state, this adiabatic 
model is also applicable to bound states. Figure [TJ compares 
the field-free potential, V g (R) +J{J+ l)/(2/iR 2 ) for J = 2, 
with the effective potentials, V g (R) +Ej(R;I), for two intensi- 
ties of the non-resonant field for 87 Rb2: as the field intensity 
increases, the height of the rotational barrier decreases. 

Equation (H)) represents an adiabatic approximation to the 
coupled rovibrational motion where the rotational wavefunc- 
tion depends parametrically on the radial variable. This ap- 
proximation is valid if the energy scales, or, respectively, 
timescales, associated with the rotational and vibrational mo- 
tion are well separated. For very large field intensities, the 
adiabatic approximation is expected to break down, because 
the matter-field interaction becomes comparable to the vibra- 
tional energy scale. This allows for coupling between states in 
different vibrational bands which is neglected in Eq. (0. 

2.3 Envisaged scheme for enhancement of photoassocia- 
tion 

Photoassociation in a magneto-optical trap can easily be sat- 
urated. Its efficiency is limited by the pair density at or near 
the Condon radius which is the internuclear distance where 
the photoassociation laser induces a resonant transition from 
the initial pair of colliding atoms to a weakly bound level 
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Fig. 1 Effective potential, 



-Ej(R;I), that the atom pair 



experiences in the presence of non-resonant laser light. The 
field-free rotational quantum number is J = 2 for which a shape 
resonance is observed at about 290 flK. The non-resonant light shifts 
the position of the shape resonance to lower energies as indicated by 
the vertical lines (resonance positions calculated with the full 2D 
Hamiltonian). 



of the atoms' electronically excited state. The ^-dependence 
of the pair density is given by Pt.i{R) = / dcos 9Pt,i{Ri 9) 
and illustrated in Fig.|2]for 87 Rb atoms. The constant behav- 
ior at large distances results from summing over many plane 
waves with random phases and reflects the equiprobability of 
finding two atoms at a certain distance. At short range, the 
interaction potential modulates the pair density. If only s- 
waves contribute significantly to the pair density, all scattering 
wave functions have their nodes at the same position at short 
range, giving rise to zeros in the pair density, cf. Pt,i{R) f° r 
T = 100/iK in Fig. [2 This reflects a threshold law— govern- 
ing collisions at very low energy. Photoassociation occurs at 
intermediate distances, typically between 40 ao and 200 ao de- 
pending on the excited state potential curve. Such a choice of 
Condon radius is a compromise between large free-to-bound 
transition probabilities (demanding large R) and sufficiently 
large binding energies required to avoid dissociation back into 
two atoms upon return to the electronic ground state (demand- 
ing short R). If formation of molecules in their electronic 
ground state proceeds via a sequence of short, optimally cho- 
sen pump and dump pulses, integrated pulse energies on the 
order of a few nano-Joule are sufficient to completely deplete 
the pair density near the Condon radius.-- Nevertheless, this 
creates only one to ten molecules per pulse sequence, depend- 
ing on the density of the atoms in the trap.— 

A further increase of the number of photoassociated 
molecules requires manipulation of the initial pair density, 
Pt,i(R)- To this end, we pursue the following idea for con- 
trolling a shape resonance with non-resonant light: A shape 
resonance is a quasi-bound state where the probability den- 
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Fig. 2 Thermal pair density p T ,t(R) = J dcos 6 p T .i{R, 0) for 87 Rb 
atoms (without any non-resonant field): at long range, the 
probability of finding two atoms at a certain distance is equally 
distributed; at short range the interaction potential leads to a 
modulation of the pair density. For T = 300 fiK, the presence of the 
close-by rf-wave shape resonance is reflected in an enhanced pair 
density at short range compared to T = 100/lK and in the nodes of 
Pt.i(R) disappearing, for example near R = 58 ao, due to states with 
J >0. 



sity of an atom pair becomes trapped inside the rotational bar- 
rier. For typical detunings, this happens at distances that cor- 
respond to the Condon radius for photoassociation 2£ and can 
therefore enhance the pair density that is available for excita- 
tion by a photoassociation laser. However, due to the thermal 
weight, e~ E " J l kBl , a shape resonance contributes significantly 
to the thermal ensemble described by Eq. (ff) only if its energy 
is close to the thermal energy of the scattering atoms. How- 
ever, the energy of shape resonances is usually of the order of 
milli-Kelvin, i.e., one to two orders of magnitude above typi- 
cal MOT temperatures. We therefore apply non-resonant light 
with an intensity chosen such that the position of the shape 
resonance is moved toward the mean energy of the scattering 
atoms. One option consists in slowly switching on the non- 
resonant laser light, thereby ensuring adiabatic following of 
the thermal cloud of atoms. The field-free eigenstates present 
in the thermal density, Eq. (0, are then transformed adia- 
batically into field-dressed states, i.e., into the eigenstates of 
Hamiltonian This modified density matrix, made up of the 
non-resonant field-dressed states, constitutes the initial state 
for photoassociation with short laser pulses.— - — Assuming 
that an optimally chosen photoassociation pulse completely 
depletes the modified pair density near the Condon radius - 
as in the field-free case— - the enhancement of the number of 
photoassociated molecules is simply given by the ratio of the 
field-dressed to the field-free thermal pair densities. 

Short-pulse photoassociation under non-resonant field con- 
trol involves two timescales: (i) a slow timescale for the 
switching of the non-resonant laser light, determined by the 



requirement of adiabaticity with respect to the rotational mo- 
tion, and (ii) a short timescale for the sequence of pump and 
dump pulses, determined by the requirement of a bandwidth 
that is optimal for photoassociation. The short timescale cor- 
responds to a few picoseconds for transform-limited pulses 
(or a few tens of picoseconds for shaped pulses). The slow 
timescale is determined by the rotational periods which are 
found to be 82ns and 36.5 ns for the rotational ground state 
of the last vibrational band for rubidium and strontium, re- 
spectively. While it is rather difficult to assign a rotational 
period to regular scattering states, the quasi-bound character 
of a shape resonance implies that its rotational period is well- 
defined. For the <f -wave resonance of rubidium and the g-wave 
resonance of strontium considered below, we find rotational 
periods of 2 /is and 350ns, respectively. A crucial question to 
be answered below is whether the modified shape resonances 
live long enough to allow for adiabatic switching of the non- 
resonant laser light. 

3 Results and discussion 

We study the non-resonant field control of a shape resonance 
for 87 Rb and 88 Sr at typical MOT temperatures, between 
50jUK and 150/lK for rubidium and about 20 /iK for stron- 
tium.— For colliding 87 Rb atom pairs, a if-wave shape res- 
onance has been observed at about 290 jUK,— and for 88 Sr, 
we find a g-wave resonance at about 1 .75 mK. Below, we in- 
vestigate the influence of non-resonant laser light on the posi- 
tion and lifetime of these two shape resonances and the result- 
ing enhancement of pair density at photoassociation distances. 
The feasibility of non-resonant field control of a shape reso- 
nance is determined by (i) the intensity required to move the 
position of the resonance close to energies corresponding to 
the MOT temperature and (ii) the lifetime of the modified res- 
onance, which must be sufficiently long to enable adiabatic 
switching of the non-resonant field. The finite lifetime of the 
shape resonance is caused by tunneling through the rotational 
barrier of the quasi-bound resonance state. Following the in- 
tuitive picture suggested by Fig. [T] where the energy of the 
field-dressed shape resonance is plotted for two laser intensi- 
ties, the lifetime of the resonance is expected to increase with 
the non-resonant field intensity. 

Calculating resonance lifetimes at very low energies is 
a challenging numerical problem.— We have therefore em- 
ployed two different methods to determine resonance lifetimes 
- complex absorbing potentials 57 and the width of peaks of 
the rotational constants, B„ = (n | \n), that occur for shape 

resonances^. Complex absorbing potentials are character- 
ized by two parameters, potential strength and width. Con- 
vergence with respect to these two parameters can be verified 
by representing the eigenvalues, (E„, ~iF n /2), in the complex 
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plane.— When the resonance occurs at very low energy, the 
potential strength needs to be small and the width huge, im- 
plying a very large spatial grid. Moreover, if the resonance is 
pushed below the dissociation threshold, the lifetime becomes 
infinite. As this limit is approached, it becomes increasingly 
difficult to achieve convergence of the parameters for the com- 
plex absorbing potential. When using the second method, the 
rotational constants, B„ = (n\ ^ 3 |n), are plotted versus en- 
ergy, E„, and fitted to a Lorentzian. 22 As the resonance posi- 
tion is shifted to lower energy by the non-resonant field, fewer 
box-discretized continuum states contribute to the peak such 
that identification of the peak becomes increasingly difficult. 
For all field intensities, the results agree with each other at 
least within 5%. This accuracy is sufficient for the order-of- 
magnitude estimate that is needed to determine the feasibility 
of the control scheme. 

3.1 Rubidium 

The potential energy curve for the lowest triplet state, a 3 £+, 
of 87 Rb2, is obtained by smoothly connecting ab initio data 
at short-range— with the asymptotic expansion C^/R 6 + 
C & /R s + C l0 /R w at long range where the C, coefficients are 
taken from Ref. [59l Polarizabilities based on ab initio calcu- 
lations are found in Ref. and fit Silberstein's formula very 
well at distances larger than 1 6 ao Q Not surprisingly, the short- 
range behavior of the polarizabilities does not play any role 
in our study, i.e., calculations with the polarizabilities from 
Ref. and calculations with polarizabilities according to Sil- 
berstein's expansion, Eq. (0, with some cutoff value at short 
distance give identical results. 

Resonance energy and lifetime as a function of the non- 
resonant field intensity / are shown in Fig. [3] for the d-wave 
shape resonance of 87 Rb, comparing the results of the 2D de- 
scription obtained by diagonalizing the Hamiltonian © with 
those obtained for the adiabatic model, Eq. (|8}. The position 
of the resonance is moved to smaller energies by increasing 
the field intensity, as expected. Overall, only moderate field 
intensities, smaller than 10 9 W/cm 2 , are required to move the 
position of the resonance close to the energy that corresponds 
to a typical MOT temperature of 100 fiK. The adiabatic ap- 
proximation underestimates the energy shift significantly for 
intensities larger than 3 x 10 8 W/cm 2 . We attribute this dis- 
agreement between the adiabatic approximation and the 2D 
description to the fact that states belonging to different vi- 
brational manifolds are mixed by a sufficiently strong non- 
resonant field. Such a mixing is not accounted for in the adia- 
batic approximation, cf. Section |2l2l 

While the adiabatic approximation, Eq. (0 captures the 
qualitative behavior of the resonance energy as a function of 



* For R > 16 ao, the relative error is below 1% (3%) for the perpendicular (par- 
allel) component. 
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Fig. 3 Energy and lifetime of the field-dressed shape resonance vs 
intensity of the non-resonant field for 87 Rb (field-free t/-wave shape 
resonance in the lowest triplet potential). Results for calculations 
based on the full 2D Hamiltonian, Eq. 10, and the adiabatic 
approximation, Eq. ([8}, are compared. The lifetimes are calculated 
employing complex absorbing potentials (CAP, ID and 2D) and 
determining the peak width of the rotational constants (B„). 

non-resonant field intensity correctly, its disagreement with 
the full 2D description, Eq. (g), with respect to the lifetimes is 
striking for the large intensities shown in Fig. [3] In fact, such 
a dramatic breakdown of the adiabatic approximation is not 
necessarily expected since overall the intensities of Fig. [3] are 
moderate and one would expect a quantitative disagreement 
such as that found for the resonance energy. The adiabatic 
approximation essentially assumes that the effect of hybridiza- 
tion due to the non-resonant field is fully captured by the mod- 
ified barrier height of effective potential, cf. Fig. [T] and it 
neglects in particular a coupling between rotational and vibra- 
tional motion. Vibrational motion that qualitatively differs for 
different rotational states is thus not correctly accounted for. 

The degree of hybridization in the resonance wave function 
is examined in Fig. @]by plotting the absolute square of the 
rotational weights, 

c'j = L[ dR [ dcosd}ife(R,e)%j(R,0), 

for the lowest four even partial waves as a function of non- 
resonant field intensity (note that the coupling mixes only par- 
tial waves of the same parity). Here, (p n j(R,9) denote the 
field-free eigenstates, i.e., the field-free box-quantized scat- 
tering states, and \ir]. es (R, 9) is the resonance wavefunction ob- 
tained by diagonalizing Hamiltonian for a given value of 
the non-resonant field intensity I. For 7 = 0, the resonance is 
a pure J = 2 state. As the non-resonant field intensity is in- 
creased, a substantial amount of / = is mixed in. Higher 
partial waves do not contribute to the resonance wavefunc- 
tion (note the different scales of the 7 = 4 and J = 6 pan- 
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Fig. 4 Contribution of different partial waves to the resonance 
wavefunction vs non-resonant field intensity for 87 Rb (field-free 
J 
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Fig. 5 Modification of the rubidium pair density: Enhancement near 
R = 75 ao (left) and ^-dependence of the field-dressed pair density, 
Lnje~ Enj/kBT J dcos6\(p nJ (R, 9)\ 2 /47tR 2 , cf. Eq. ©, for zero field 
(grey), / = 6 x 10 9 W/cm 2 (green) and / = 9 x 10 9 W/cm 2 (blue). 



els). A turnover in the resonance lifetime as a function of 
field intensity is observed in the right-hand side of Fig. [3] At 
the corresponding intensity, the 7 = contribution amounts to 
about 30%, and for the largest intensity shown in Fig. |4] the 
resonance wavefunction is predominantly of 7 = character. 
The different behavior of the lifetime in the full 2D Hamilto- 
nian and the adiabatic approximation is now easily rational- 
ized by the qualitatively different /^-dependences of the 7 = 2 
and 7 = states: While the field-free 7 = 2 shape resonance 
is a quasi-bound state confined to short range by the centrifu- 
gal barrier, all 7 = states are of a purely scattering character. 
As the non-resonant field intensity is switched on, the reso- 
nance wavefunction becomes a superposition of these states, 
effectively loosing its quasi-bound character. 

The effect of the non-resonant field on the pair density is 



shown in Fig. [5] The pair density at photoassociation dis- 
tances is enhanced by about one order of magnitude, cf. left- 
hand side of Fig. [5] due to increasing the thermal weight of 
the shape resonance by moving its energy close to kg times 
the trap temperature. Nevertheless, the enhancement is larger 
for T = 50/iK compared to T = 100 /xK although the reso- 
nance position is only moved to T = 100 /J.K, cf. Fig. [3] for 
the largest intensity shown in Fig. This is due to the fact 
that for T = 100/iK, already the field-free pair density is in- 
fluenced by the presence of the shape resonance, while its ef- 
fect on the field-free pair density is negligible for 7 = 50/iK. 
The change brought about by the non-resonant field is thus 
larger for T = 50/xK. The enhancement of the pair density 
at all short distances is illustrated by the right-hand side of 
Fig. [5] which shows the /^-dependence of the field-dressed 
pair density for three different intensities of the non-resonant 
field. The plotted quantity corresponds to the unnormalized 
pair density where the 9 -dependence has been integrated over. 
While a few partial waves are sufficient to converge the R- 
dependence of the pair density, the normalization factor re- 
quires a much larger number of partial waves for convergence . 

The pair density is significantly enhanced also at distances 
shorter than the last peak of the last bound level. However, 
due to the fast oscillations of the wavefunction, this pair den- 
sity enhancement does not translate into substantially larger 
number of molecules that can be photoassociated. The largest 
increase of the photoassociation efficiency is due to the pair 
density enhancement close to the last peak of the last bound 
level, i.e., about 75 ao for rubidium. 

In summary, applying a non-resonant laser field shifts the 
position of the of-wave rubidium shape resonance to smaller 
energies. While this leads to an increasing lifetime of the res- 
onance in an effective ID model, the shape resonance looses 
its quasi-bound character in the full 2D description since it is 
mixed with states of pure scattering character. The lifetime 
therefore decreases after an initial increase. The rubidium 
atom pair density at distances relevant to photoassociation is 
found to increase by about one order of magnitude. 

3.2 Strontium 

The potential energy curve for the electronic ground state of 
88 Sr2 has been obtained spectroscopically, and an analytical fit 
was reported in Ref. |6j] We approximate the polarizabilities 
by their long range expansion, Eq. (0, with the atomic values 
taken from Ref. |62| 

Resonance energy and lifetime as a function of the nonres- 
onant field intensity / are displayed in Fig. [6] for the g-wave 
shape resonance of 88 Sr, comparing again results obtained for 
the 2D Hamiltonian (0]i with those of the adiabatic approxi- 
mation, Eq. ([S]). Like in the case of rubidium, cf. Fig. [3] the 
position of the resonance is moved to smaller energies by in- 
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non-resonant field intensity ( 10 W/cm ) 

Fig. 6 Energy and lifetime of the field-dressed shape resonance vs 
intensity of the non-resonant field for 88 Sr (filled symbols: field-free 
g-wave (7 = 4) shape resonance in the electronic ground state 
potential, empty symbols: field-free 7 = 8 resonance). Results for 
calculations based on the full 2D Hamiltonian, Eq. {4])> and the 
adiabatic approximation, Eq. ([§}, are compared. The lifetimes are 
calculated employing complex absorbing potentials (CAP, ID) and 
determining the peak width of the rotational constants (B„, ID and 
2D). Note the logarithmic scale for the lifetimes. 



creasing the non-resonant field intensity. Slightly larger field 
intensities are required than for rubidium. This is easily ra- 
tionalized in terms of the different atomic polarizabilities and 
rotational constants. The atomic polarizability of rubidium is 
larger than that of strontium, yielding a stronger interaction. 
Moreover, the rotational constant of strontium is larger than 
that of rubidium, so in order to achieve the same effect, a larger 
non-resonant field is required. Nevertheless, the intensities re- 
main moderate for strontium as well: about 5 x 10 9 W/cm 2 
are required to move the position of the resonance close to 
the energy corresponding to a MOT temperature of 20 fiK, 
typical for the two-color MOTs employed for alkaline-earth 
species.-^, Also, like for rubidium, the adiabatic approxima- 
tion significantly underestimates the energy shift at the large 
laser intensities. This becomes particularly evident at inten- 
sities in excess of 3 x 10 9 W/cm 2 . Inspecting the pair density 
enhancement for strontium (see Fig. |8]below), a second res- 
onant feature is observed which is attributed to a field-free 
7 = 8 shape resonance. The position of this resonance vs 
non-resonant field intensity is traced by empty red circles in 
Fig. [6] While the field-free 7 = 8 resonance starts out at an 
energy much larger than that of the 7 = 4 resonance, its en- 
ergy is more strongly affected by the non-resonant field than 
that of the 7 = 4 one. It is therefore moved to lower tem- 
peratures faster. The positions of the two resonances cross at 
about 3.7 x 10 9 W/cm 2 . Since both states have the same sym- 
metry, the shape resonances exhibit an avoided crossing as I 
is varied. This field-induced phenomenon is characterized by 
a strong coupling between the vibrational and rotational mo- 
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Fig. 7 Contribution of different partial waves to the resonance 
wavefunction vs non-resonant field intensity for 88 Sr (field-free 
7 = 4 shape resonance on the left-hand side, field-free 7 = 8 shape 
resonance on the right-hand side). 

tion.— It cannot be described by the adiabatic model which 
treats these two degrees of freedom as independent. 

While for both rubidium and strontium the field-free life- 
times are on the order of 50 ns, a dramatic increase of the 
lifetime is observed for strontium as the non-resonant field 
is applied (note the logarithmic scale of the lifetimes on the 
right-hand ordinate in Fig. |6). This is in a stark contrast with 
rubidium, cf. Fig. [3] where an increasing lifetime is observed 
only in the adiabatic approximation, but was not confirmed by 
the 2D calculation. In strontium, the 2D lifetimes decrease 
with the non-resonant field intensity in an intermediate inten- 
sity range, but as the intensity is further increased, the life- 
times go up by an order of magnitude. The dip in the life- 
times, Fig. [6] occurs in a range of non-resonant field intensi- 
ties where the 7 = 4 and 7 = 8 resonances are close in energy. 
Past the avoided crossing, the lifetime of the field-free 7 = 4 
resonance increases again due to the mixing that takes place 
in the crossing region: The 7 = 8 resonance is quite narrow 
and thus long-lived. Due to the interaction of the two reso- 
nances, this character of the resonance is partially transferred 
to the field-free 7 = 4 resonance. Despite the fact that both 
the adiabatic approximation and the 2D description predict an 
overall increase of the lifetimes, the adiabatic approximation 
ceases to be valid for strontium (like for rubidium) at interme- 
diate and large intensities since quantitatively both energies 
and lifetimes do not agree with those of the 2D description: 
the more conservative estimates of the lifetime obtained in the 
adiabatic approximation by determining the peak width of the 
rotational constants are still an order of magnitude larger than 
the 2D values. As for rubidium, the breakdown of the adia- 
batic approximation is attributed to mixing of states from dif- 
ferent vibrational manifolds. 

A mixing of states from different vibrational manifolds 
leads to particularly pronounced effects if two resonances 
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Fig. 8 Modification of the strontium pair density: Enhancement 
near R = 55 ao (left) and ^-dependence of the field-dressed pair 
density, £„ y e - E ><j/ k i> T Jdcos 6\<p nJ {R, 8)\ 2 /4xR 2 , cf. Eq. l©, for 
T = 20,uK (right). The lightblue and grey symbols in the left panel 
display the averaged enhancement (pair density integrated from 
R = 32 ao to R = 120ao.) Note the different scales of the right-hand 
side panels. 



are concerned, i.e., for non-resonant field intensities close to 
1 = 3.1 x 10 9 W/cm 2 , cf. left-hand side of Fig. This is il- 
lustrated in Fig. [7] in terms of the partial contributions to the 
resonance wavefunctions for both the field-free 7 = 4 reso- 
nance (left-hand side of Fig. [7]) and the field-free 7 = 8 shape 
resonance (right-hand side of Fig.|7J. For 7 = 0, the resonance 
wavefunctions are pure 7 = 4 and 7 = 8 states, respectively. 
As the non-resonant field is switched on, a substantial amount 
of first 7 = 2 and then 7 = is mixed in to the field-free 7 = 4 
state. While this resonance wavefunction is predominantly of 
7 = 2 character at intermediate intensities, it acquires essen- 
tially an s- wave character at large intensities. Close to the in- 
tensity where the two resonances cross, some 7 = 8 character 
is mixed in as well (little orange peak in the left-hand side of 
Fig. [7}. The presence of the second resonance shows up in all 
the partial wave contributions causing discontinuous behavior 
near the crossing point. Partial waves corresponding to 7 = 10 
and higher do not contribute to the resonance wavefunctions 
(grey symbols in the left-hand side of Fig. |7). The field-free 
7 = 8 shape resonance (right-hand side of Fig. [7} essentially 
keeps its character with just a small admixture of 7 = 6 until 
the crossing region with the 7 = 4 resonance is reached. There 
it acquires a substantial amount of 7 = 2 and 7 = character. 
This is due to the strong mixing between both states induced 
by the avoided crossing. Once the crossing region is passed, 
the resonance wavefunction recovers its mainly 7 = 8 charac- 
ter. This represents another unambigious signature of interac- 
tion between the two resonances. 

The presence of the two resonances shows up clearly also 
in the modification of the pair density, cf. Fig. [8] While the 



first peak in the enhancement (left-hand side of Fig. |8), near 
I = 4 x 10 9 W/cm 2 , is due to the field-free 7 = 8 resonance 
whose energy approaches that corresponding to the trap tem- 
perature, the second peak, near I = 5 x 10 9 W/cm 2 , can be 
attributed to the field-free 7 = 4 resonance. Unlike rubidium, 
where the nodal structure of the pair density was only slightly 
modified by shifting down the resonance, cf. Fig. [5] the R- 
dependence of the pair density depends sensitively on the non- 
resonant field intensity, cf. right-hand side of Fig. [8] The en- 
hancement at a specific position is therefore compared to an 
averaged enhancement obtained by integrating the pair den- 
sity in all relevant photoassociation distances, from R = 32 ao 
to R = 120 ao (light-blue and grey symbols in the left-hand side 
of Fig. [8}. The overall trend and order of magnitude behavior 
of the averaged enhancement is the same as the enhancement 
at a specific distance, R = 55 ao. The sensitive /^-dependence 
of the pair density on laser intensity reflects the contribution 
of the two resonances which, in an effective description, corre- 
spond to two different centrifugal barriers with different turn- 
ing points. It is the outer turning point at each of the cen- 
trifugal barriers that causes the large peaks near R = 45 ao for 
/ = 4x 10 9 W/cm 2 and near R = 60 a for / = 5 x 10 9 W/cm 2 . 
Note that the pair density is also enhanced at the laser inten- 
sities where the two resonances exhibit the avoided crossing, 
/ = 3.75 x 10 9 W/cm 2 . However, the effect is most visible for 
a temperature around T = 200 jiK, which is much larger than 
typical strontium MOT temperatures. 

In summary, the interaction of two resonances leaves a pro- 
nounced signature in the atom pair density, one that should be 
observable in an experiment which probes the pair density.™ 
Moreover, for zero-field resonances which are far from the 
trap temperature, a large enhancement of the atom pair den- 
sity - here about 2 orders of magnitude - and an increase in 
the field-free resonance lifetime are found. 

We note that an enhanced pair density at photoassociation 
internuclear separations brought about by tuning a shape res- 
onance could also enable photoassociation based on Stimu- 
lated Raman Adiabatic Passage (STIRAP) 65 either within a 
single pair of pump and Stokes pulses or within a sequence of 
phase-locked STIRAP pulse pairs.^i The feasibility of STI- 
RAP photoassociation depends on a sufficient isolation of the 
initial state from the scattering continuum. A possibility to 
achieve this discussed in the literature is based on utilizing a 
Feshbach resonance.— If, however, no resonance is present, 
i.e., for an unstructured continuum, STIRAP fails. There- 
fore, experiments that made use of STIRAP for transferring 
molecules to their ground state started from molecules that 
had been already associated via a Feshbach resonance. 1 i 2 i 69 i 7Q 
A shape resonance that is brought to the right energy by non- 
resonant laser light might provide an alternative way to isolate 
the initial state for STIRAP from the scattering continuum. It 
would not rely on the presence of Feshbach resonances which 
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do not occur for example for the even isotope species of alka- 
line earth atoms. 7 However, STIRAP induces a second slow 
timescale, due to the requirement of adiabaticity with respect 
to the initial state. This translates into an even longer lifetime 
of the shape resonance that is required to guarantee success 
of STIRAP photoassociation compared with the non-adiabatic 
picosecond pump-dump photoassociation sequence discussed 
above. 

4 Conclusions 

We have considered the interaction of non-resonant laser light 
with pairs of colliding atoms held in a magneto-optical trap, 
employing a rigorous treatment of the thermal ensemble of 
atoms and comparing a 2D description of the coupled rovi- 
brational motion to an effective ID model based on the adia- 
batic approximation. Specifically, we have studied the influ- 
ence of the non-resonant light on the position and the lifetime 
of shape resonances. Such quasi-bound states, that are trapped 
behind the centrifugal barrier imposed by the angular momen- 
tum involved in the collision, lead to significant enhancement 
of the pair density at comparatively short range. This en- 
hancement readily translates into an increased efficiency of 
molecule-formation processes such as photoassociation. 20 ' 21 
Since the photoassociation rate is limited by the pair density 
at or near the Condon radius, 38,53 utilizing shape resonances 
could overcome the main obstacle toward forming large num- 
bers of molecules. 

However, most of the time, the thermal weight of shape res- 
onances at MOT temperatures is very small because their en- 
ergy is too high. In the present study, we find that the energy 
of shape resonances can be decreased over several orders of 
magnitude by applying a non-resonant laser field. It is the 
second order nature of the laser-matter interaction that guar- 
antees a monotonous decrease in energy. Moderate intensi- 
ties of the non-resonant field are sufficient: intensities of less 
than 10 9 W/cm 2 are required to move the d-wave resonance 
of rubidium from 300 jJ.K to below 100 fiK, and intensities of 
about 5 x 10 9 W/cm 2 shift the g-wave resonance of strontium 
from 1 .75 mK to below 20 jiK. The value of the required non- 
resonant field intensity is determined by the rotational con- 
stant and the atomic polarizability of a given species: small 
rotational constants and large polarizabilities are desirable to 
maximize the matter-field interaction and minimize the (com- 
peting) rotational kinetic energy. 

As the non-resonant field intensity is increased, the lifetime 
of a shape resonance first increases, then drops again. We ra- 
tionalize the increase within the ID picture by the increase in 
the tunneling time through the centrifugal barrier which ac- 
companies the lowering of the energy of the resonance. How- 
ever, the ID picture captures only part of the story: as the field 
intensity is further increased, strong hybridization occurs. As 



scattering states belonging to different partial waves get mixed 
in, the resonance wavefunction loses its quasi-bound charac- 
ter, whereby tunneling out of the centrifugal barrier is actually 
enhanced. This is the generic behavior of an isolated shape 
resonance which is observed in our example of the J = 2 res- 
onance in rubidium. A 2D description of the coupled rovibra- 
tional motion is essential to capturing this effect correctly even 
for comparatively moderate intensities. The adiabatic approx- 
imation, which assumes the vibrational and rotational degrees 
of freedom to be independent for all states, ceases to be valid. 

The physical situation changes once more if more than one 
resonance comes into play, as in our example of strontium 
where a broad J — A resonance and a narrow 7 = 8 resonance 
interact. For strong fields, the two resonances show an avoided 
crossing, in the course of which part of the narrow -resonance 
character of the 7 = 8 resonance is transferred to the J = A 
resonance, increasing the lifetime of the latter. Of course, the 
interaction of different resonances also requires a 2D descrip- 
tion of the coupled rovibrational motion and is not accounted 
for in an effective ID treatment. 

For both single and interacting resonances, significant en- 
hancement of the thermal pair density at short distances is 
found. The magnitude of the enhancement depends on the dis- 
tance between the field-free resonance position and the trap 
temperature. Maximum enhancement is obtained when the 
non-resonant field intensity is chosen such that the position 
of the shape resonance matches the trap temperature. In our 
examples, an enhancement on the order of 10 and 100 was 
found for rubidium and strontium, respectively. Such an en- 
hancement can be probed directly by pulses that are short on 
the timescale of the rotational and translational motion in the 
trap, for example by picosecond pulses. 64 

Our thermally averaged pair densities describe the struc- 
ture of the rovibrational eigenvalue problem for a given non- 
resonant field intensity. This gives rise to both a dynamic 
and a static regime that can be pursued to implement non- 
resonant field control of shape resonances in an experiment. 
The dynamic regime consists of applying a non-resonant laser 
pulse which is long compared with the rotational motion as- 
sociated with the shape resonance. The shape resonance adi- 
abatically follows the pulse and is shifted toward lower ener- 
gies. Once the energy corresponding to the trap temperature 
is reached, a resonant photoassociation or probe pulse is ap- 
plied that is short with respect to the non-resonant laser pulse. 
This scheme is expected to be successful if the timescale of 
the adiabatic switching is of the order of the lifetime of the 
resonance or shorter. In our study, we found rotational periods 
of 2 /is and lifetimes of up to 100ns for rubidium, whereas 
for strontium we found rotational periods of 350 ns and life- 
times of 500ns. Therefore, adiabatic following cannot be en- 
sured for rubidium. While losses due to the lifetime of the 
resonance might hamper the realization of a full two-orders - 
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of-magnitude enhancement, the enhancement found for stron- 
tium is significant. In general, controlling a shape resonance 
with non-resonant laser light based on adiabatic following will 
be successful for species with short rotational periods (or large 
rotational constants). 

As an alternative to adiabatically shifting the position of a 
shape resonance, a static scheme can be pursued. This is based 
on the fact that the thermally averaged atom pair density de- 
scribes a thermal equilibrium in the trap for a given temper- 
ature and non-resonant field intensity. The static scheme in- 
volves application of a constant non-resonant field and a hold 
period for thermaliziation of the atom cloud in the presence 
of the field. The field-modified equilibrium can then serve as 
a starting point for photoassociation, using either continuous- 
wave lasers or laser pulses. 

Finally, we wish to point out that our approach differs 
from previous studies on the resonance-enhancement of the 
short-range pair density and the subsequent molecule forma- 
tion—"— in that we seek to actively control the resonance. 
This is in contrast to the literature which discusses the in- 
fluence of a given resonance on the desired process, e.g., 
molecule formation. Following the spirit of coherent control 
where resonances can be used to imprint a well-defined phase 
onto the initial state of the desired process,— in our present 
study the resonance itself is the subject of control, in order 
to afford an optimal use of it. In terms of controllability, our 
main finding is that shape resonances can be controlled in their 
energetic position but not at the same time in their vibrational 
character (purely scattering vs quasi-bound). This insight de- 
fines the challenge for future work - to identify a fully con- 
trollable scattering resonance. 
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